Final  Report  for  AOARD  Grant  FA23 86-09- 1-4027 


Homeland  security  enforcement  using  novel  terahertz  technology 
December  1,  2010 

Name  of  Principal  Investigators:  Adrian  Dobroiu 
e-mail  address  :  dobroiu@riken.jp 
Institution  :  RIKEN 
Mailing  Address  : 

RIKEN,  Terahertz  Sensing  and  Imaging  Laboratory 
519-1399  Aoba,  Aramaki,  Aoba-ku,  Sendai  980-0845,  Japan 

-  Phone:+81-22-228-2124 

-  Fax:+81-22-228-2128 

Period  of  Performance :  06/03/2009  -  1 0/02/20 1 0 


Abstract:  The  purpose  of  this  research  was  to  investigate  the  application  an  ellipsometric  technique  in 
the  determination  of  the  optical  properties  of  liquids  in  the  terahertz  frequency  range.  The  main  results 
obtained  in  this  period  were  as  follows.  1.  The  reasons  for  the  previous  problems  with  repeatability 
were  identified  and  the  problems  were  largely  solved.  They  consist  in  the  occurrence  of  a  Fabry-Perot 
effect  within  the  optical  setup,  which  makes  the  signals  sensitive  to  very  slight  misalignments.  2. 
Other  phenomena  that  could  affect  the  results  by  changing  the  polarization  state  of  the  radiation  were 
found,  explained,  and  canceled.  They  relate  partly  to  the  same  Fabry-Perot  effects  and  partly  to  the 
changes  in  the  polarization  state  due  to  the  passage  of  a  focused  beam  through  circular  apertures 
whose  diameter  is  comparable  to  the  wavelength.  3.  Advances  were  made  in  the  signal  and  data 
processing.  Instead  of  complex  analytical  calculations  we  chose  to  realize  a  computer  program  that 
simulates  the  beam  propagation  and  interaction  with  the  various  elements  of  the  setup.  This  allows  a 
fitting  procedure  to  find  the  values  of  the  optical  properties  of  the  sample  liquid  by  adjusting  a  set  of 
simulation  parameters  directly  related  to  physical  parameters  in  the  actual  setup.  4.  Investigations  at 
other  terahertz  frequencies  showed  that  such  measurements  have  the  potential  to  bring  in  more  data 
which  would  allow  an  even  more  reliable  identification  of  the  sample  liquid,  but  that  for  such 
measurements  to  contribute  significantly  it  is  necessary  to  use  terahertz  frequencies  more  widely 
spaced  than  available  from  the  BWO  source  that  we  used  in  our  experiments.  5.  By  using  a  larger 
silicon  prism,  the  errors  were  greatly  reduced  as  the  beam  was  used  more  efficiently. 

Introduction:  Terahertz  radiation  is  the  electromagnetic  waves  with  a  frequency  around  1  THz  = 
10 12  Hz.  In  the  electromagnetic  spectrum,  this  range  is  found  between  high  frequency  millimeter 
waves  and  long  wavelength  infrared  radiation.  Until  about  15  years  ago  the  terahertz  frequency  range 
was  very  little  studied  and  material  properties  in  this  range  were  largely  unknown,  because  there  were 
not  many  affordable  sources  and  detectors  to  work  at  these  frequencies.  As  soon  as  these  became 
available,  the  research  in  this  field  exploded,  so  that  now  there  are  hundreds  of  scientific  papers 
published  every  year  just  concerning  this  frequency  range  and  many  practical  applications  have 
appeared  which  prove  that  it  has  a  very  good  potential  to  help  solve  real-world  problems:  post-office 
screening  of  envelopes  to  detect  illicit  drugs,  airport  checks  for  detection  of  weapons  and  explosives 
hidden  under  passenger  clothes,  diagnosis  of  cancer,  bums  and  other  medical  conditions,  art 
conservation  and  diagnostics,  agricultural  automation,  industrial  investigation,  etc.  The  main  qualities 
that  recommend  this  frequency  range  are  the  fact  that  many  otherwise  opaque  materials  are 
transparent  to  terahertz  waves,  the  fact  that  many  organic  materials  have  specific  fingerprint  spectra 
allowing  their  identification,  and  the  fact  that,  unlike  for  instance  X  rays,  terahertz  waves  have  no 
detectable  impact  on  the  living  tissue,  which  makes  them  safe  to  use  where  operator  health  is  a 
concern. 

The  ultimate  goal  of  our  research  is  to  develop  a  technique  for  identifying  liquids.  This  can  have 
numerous  applications,  from  homeland  security  (for  instance  checks  in  airports)  to  pharmaceutical  and 
chemical  analysis. 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

05  OCT  2010 


2.  REPORT  TYPE 

Final 


3.  DATES  COVERED 

03-06-2009  to  02-07-2010 


4.  TITLE  AND  SUBTITLE 

Homeland  Security  Enforcement  Using  Novel  Terahertz  Technology  V 


5a.  CONTRACT  NUMBER 

FA23860914027 


6.  AUTHOR(S) 

Adrian  Dobroiu 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS (ES) 

RIKEN  519-1399, Aoba,  Aramaki,  Aoba-ku, Sendai 
980-0845, Japan, NA,NA 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS  (ES) 

AOARD,  UNIT  45002,  APO,  AP,  96338-5002 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AOARD 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AOARD-094027 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

This  is  a  report  of  research  on  the  application  of  an  ellipsometric  technique  to  determine  the  optical 
properties  of  liquids  in  the  terahertz  frequency  range. 

15.  SUBJECT  TERMS 

Terahertz  Technology,  Plastics  Explosives,  Biological  Agents 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

6 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


The  principle  of  liquid  identifying  is  to  allow  the  interaction  between  a  terahertz  beam  and  the 
unknown  liquid  to  interact,  which  gives  information  about  the  absorption  coefficient  and  refractive 
index  of  the  liquid.  These  two  optical  properties  are  different  from  one  liquid  to  another. 

If,  by  an  unlikely  coincidence,  two  liquids  have  the  same  pair  of  optical  properties  at  one  frequency, 
then  they  can  most  likely  be  discriminated  at  another  frequency.  We  deemed  it  to  be  practically 
impossible  to  create,  even  intentionally,  a  mixture  of  liquids  that  imitate  the  optical  properties  of  a 
given  liquid  at  all  terahertz  frequencies.  As  a  consequence  we  believe  that  a  machine  employing 
terahertz  radiation  to  identify  liquids  cannot  be  fooled. 

This  means  that,  by  using  a  database  of  liquid  optical  properties,  measuring  the  absorption  coefficient 
and  the  refractive  index  of  an  unknown  liquid  will  allow  the  reliable  identification  of  the  unknown 
liquid. 

Experiment:  The  experimental  setup  consists  of  several  basic  elements:  a  terahertz  source  and  an 
appropriate  detector,  a  silicon  prism  which  allows  the  interaction  between  the  radiation  and  the  liquid 
in  a  controlled  way,  and  a  pair  of  polarizers,  as  shown  in  Figure  1 . 


Sample 


Figure  1 .  Schematic  of  the  experimental  setup,  showing  only  the  main  components.  The  optical  chopper  placed  in 
front  of  the  BWO  source  is  non-essential;  it  is  used  to  achieve  a  low-noise  detection. 


The  first  polarizer  is  fixed  so  that  it  makes  45°  with  both  axes,  while  the  second  polarizer  (labeled 
“analyzer”)  is  mounted  on  a  computer-controlled  motorized  rotation  stage  and  is  used  to  analyze  the 
polarization  properties  of  the  radiation  after  it  interacted  with  the  liquid  sample.  The  same  computer 
that  controls  the  rotation  stage  reads  and  processes  the  signal  from  the  detector  synchronously  with  the 
analyzer  rotation. 

The  prism  used  in  the  setup  is  made  of  high-resistivity  silicon,  which  is  the  most  transparent  material 
available  for  the  terahertz  range.  Compared  to  the  prism  we  used  in  previous  measurements,  this  one 
is  large  (the  small  sides  are  about  55  mm)  and  allows  an  efficient  use  of  the  available  radiation,  with 
minimum  losses.  The  upper  surface  of  the  prism  has  an  edge  rail  made  of  silicone  rubber  that  prevents 
the  liquid  from  spilling  over  the  edge.  A  correct  interaction  between  the  terahertz  radiation  and  the 
liquid  requires  a  liquid  thickness  that  is  considerably  larger  than  the  so-called  penetration  depth  of  the 
liquid.  This  penetration  depth  differs  from  liquid  to  liquid,  but  is  generally  in  the  order  of  tenths  of  a 
millimeter;  consequently  we  chose  a  liquid  thickness  of  about  3  mm. 


We  made  experiments  on  several  liquids:  water,  alcohol,  various  water-alcohol  mixtures,  acetone, 
milk,  etc.  The  typical  curves  that  were  obtained  are  shown  in  Figure  2. 
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Figure  2.  Examples  of  signals  obtained  on  various  liquids. 
The  water-ethanol  mixture  is  50%-50%  by  mass. 


The  graph  shows  that  at  least  for  these  sample  liquids  there  is  no  overlapping  and  the  liquids  can  be 
easily  identified  even  without  further  processing  of  the  data. 

For  easy  understanding,  a  2D  map  is  shown  in  Figure  3,  where  each  liquid  is  represented  as  a  dot. 
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Figure  3.  A  2D  representation  of  the  data.  Air  data  is  also  given  for  reference.  The  liquid  labeled 
E5+W95  is  a  mixture  of  5%  ethanol  in  water  by  mass  (a  concentration  similar  to  that  of  beer). 
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The  optical  setup  contains  elements  used  for  alignment,  such  as  a  visible  light  source  that  creates  a 
beam  exactly  overlapping  the  terahertz  beam,  a  removable  small  aperture  used  for  checking  if  the 
prism  has  moved.  Occasionally  we  used  a  partial  absorber  (attenuator)  to  cancel  the  Fabry-Perot 
effects  caused  by  repeated  reflection  between  the  source  and  the  detector. 

Data  processing:  The  experimental  data  were  processed  using  a  computer  program  that  we  made  for 
this  purpose.  In  essence,  the  program  simulates  the  experimental  setup,  allowing  for  the  adjustment  of 
many  parameters  such  as  the  polarizer  angle,  an  offset  of  the  analyzer  angle,  tilt  angle  for  the  silicon 
prism,  errors  of  the  prism  angles,  the  refractive  index  of  the  prism  and  that  of  the  liquid,  etc. 


By  repeatedly  adjusting  these  parameters  the  program  allows  the  fitting  of  the  measured  signal  with 
the  simulated  signal.  This  gives  the  two  optical  parameters  of  the  unknown  liquid  as  an  output  of  the 
program. 

In  principle  the  same  result  should  be  obtained  analytically,  by  solving  the  equations  governing  the 
setup.  However,  the  complexity  of  the  setup  and  the  fact  that  the  radiation  suffers  a  total  of  three 
Fresnel  interactions  (two  refractions  and  one  reflection)  makes  such  analytical  calculation  extremely 
difficult  and  any  additional  parameter  complicates  the  formula  beyond  usability.  We  decided  that  the 
simulation  option  was  the  best  available. 

Results:  The  first  result  is  the  experimental  confirmation  that  two  liquids  give  two  different  output 
signals,  which  is  the  most  important  point  in  any  identification  process.  By  further  data  processing  we 
were  able  to  obtain  the  optical  parameters  of  the  liquids.  This  allows  the  use  of  data  measured  and 
published  by  others,  so  that  we  don’t  have  to  measure  everything  by  ourselves  in  order  to  make  the 
database. 

Besides  this  main  result,  we  believe  that  several  other  observations  made  during  the  experiments  are 
interesting  both  as  theoretically  and  practically.  An  important  observation  was  the  fact  that 
Fabry-Perot  resonances  inside  the  optical  path  have  a  strong  negative  effect  on  the  measurements  and 
must  be  canceled.  They  are  particularly  destructive  in  ellipsometry  measurements  because  the  quality 
factor  of  the  resonances  changes  with  the  rotation  of  the  analyzer;  this  has  an  impact  on  the  perceived 
linearity  of  the  detector  and  eventually  on  the  measured  data.  The  effect  is  so  strong  that  the  procedure 
we  used  to  fit  the  ellipsometric  data  with  a  sine  wave  can  no  longer  work,  as  it  gives  negative  values 
for  part  of  the  fitted  data,  which  is  impossible  (any  radiation  intensity  must  be  non-negative). 

In  order  to  cancel  the  resonances  we  had  two  options.  One  is  to  use  an  attenuator  somewhere  along  the 
optical  path.  The  absorption  must  be  strong  enough  to  sufficiently  extinguish  the  beam  coming  back 
after  a  reflection,  but  weak  enough  to  let  some  radiation  pass.  The  other  option  is  to  modulate  the 
frequency  of  the  terahertz  source  so  that  the  effects  of  many  frequencies  average  in  a  way  that  the 
effects  of  the  Fabry-Perot  resonance  are  canceled.  By  far  the  best  option  is  the  latter  one,  because  it 
ensures  a  good  use  of  all  the  terahertz  power  available,  but  the  first  option  must  be  retained  for  the 
case  where  in  a  practical  implementation  a  source  is  used  that  does  not  allow  frequency  modulation. 

Another  useful  observation  was  the  fact  that  the  small  aperture  we  used  for  checking  the  prism 
alignment  had  an  effect  on  the  polarization  state  of  the  beam,  in  that  it  reduced  the  polarization 
contrast  of  a  linearly  polarized  beam.  We  still  have  to  investigate  what  the  theoretical  reason  for  that 
observation  may  be,  but  from  a  practical  point  of  view  we  found  that  it  is  important  to  remove  the 
aperture  before  the  data  is  measured.  Another  implication  of  this  is  that  we  cannot  use  a  detector 
whose  entrance  pupil  is  so  small  as  to  affect  the  polarization. 

Discussion:  The  research  presented  here  shows  that  in  principle  it  is  possible  to  identify  liquids. 
Before  this  principle  can  be  implemented  in  a  practical  application,  however,  it  is  necessary  to  solve  a 
few  problems.  Among  them,  probably  the  most  challenging  one  is  to  achieve  measurements  while  the 
liquid  is  still  in  its  container.  This  means  that  the  effect  of  the  container  on  the  signal  must  be  either 
known  or  somehow  canceled.  One  possible  solution  for  this,  that  we  haven’t  tried  out  yet,  is  to  use  a 
pulsed  terahertz  beam,  which  would  allow  a  gated  measurement,  that  is,  the  data  would  be  collected 
only  from  the  interaction  of  the  beam  with  the  liquid,  not  with  the  container. 

However,  the  research  can  be  used  in  its  current  stage  for  other  applications  such  as  the  measurement 
of  concentration  of  mixtures.  For  this  purpose  a  series  of  mixtures  with  known  concentrations  are 
measured,  then  the  place  the  mixture  with  the  unknown  concentration  is  found  among  the  others  by 
interpolation.  As  a  rough  approximation  it  is  possible  to  measure  just  the  two  pure  liquids  and  assume 
that  the  optical  properties  of  the  mixture  are  linearly  distributed  between  the  extremes,  but  from  our 
experience  with  mixtures  of  water  and  ethanol  the  error  such  an  assumption  leads  to  are  too  large  for 
any  practical  purpose,  so  at  least  a  few  intermediate  concentrations  must  be  checked. 


List  of  Publications:  For  now  only  the  work  below  was  published  on  this  subject,  presented  at  the 


annual  conference  on  terahertz  physics.  We  plan  on  publishing  a  full  journal  paper  soon. 


•  Adrian  Dobroiu  and  Chiko  Otani,  “Measurement  of  the  complex  refractive  index  of  liquids  in  the  terahertz 
range  using  ellipsometry,”  35th  International  Conference  on  Infrared,  Millimeter  and  Terahertz  Waves 
(IRMMW-THz  2010)  Sept.  5-10,  2010,  Rome,  Italy,  (paper  Tu-P.51) 

A  copy  of  the  poster  presented  at  the  conference  is  attached  on  the  next  page. 


Measurement  of  the  complex  refraction  index  of  liquids 
in  the  terahertz  range  using  ellipsometry 

Adrian  Dobroiu.  Chiko  Olani  ikmm\v-iiii2oiu 

klKI-'N  —  Terahertz  Sonxinj*  and  Im»j»!nj»  I  ahnnitnry  Sendai  Japan 


Purpose 


■  Liquid  identification 

■  liquid  anulvsis 

*  Solution  coixxntratian  measurement 

the  technique  could  is  applied  wherever  we  need  to  identify  or  analyze 
liquids,  for  example  in  homeland  security,  chemistry,  fond  production, 

phuTmtcv.  etc 


We  have  described  a  similar  technique, 
but  with  a  afferent  pnnciple.  in: 

A.  Dohrma  R.  Reigang.  C.  (Xani.  and 
k.  kawase.  "Monolithic  fabry -Perot 
resonator  fur  the  measurement  of  optical 

constants  ir  the  terahertz  range"  Applied 
Physics  Letters  86. 231 107  ( 2005 ) 


Principle 

A  terahertz  wave  is  sent  through  a  transparent  silicon  prism.  W ithout  tire  liquid 

sample  a  total  reflection  occurs  on  the  upper  surface.  When  a  liquid  sample  is 

placed  on  top  of  the  prism,  the  total  internal  reflection  becomes  partial  reflection, 
and  the  intensity  and  polanzation  state  of  the  wave  are  changed,  winch  allows 
measuring  the  the  complex  retraction  index  of  the  sample 

The  tfHRmftl  reflecTion  is  governed  hv  two  Frennel  reflection  iwttfeients  one  tor 
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•  the  m's  arc  the  complex  refraction  indices, 
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lhe  retract  rvc  index  of  the  sample  can  he  measured  ha'  analyzing  the  polarization 
suit  ol' lhe  output  beam.  We  propose  using  an  dlipsoroculc  technique 

1 .  lhe  input  beam  is  polarized  al  4>".  halfw  ay  between  the  parallel  and  the 
perpendicular  polarization. 

2.  l  he  output  beam  is  passed  through  a  second  polarizer,  w  hose  orienuiion  can 
be  adujsted. 

1  Bv  .*  uLluyj  lhe  second  pulai  i/et.  a  sine-shapeil  inlensilv  sinna!  is  teuiukd 
I  -  It  (l  4  a  ■  cos  2a + b  -  k in  2a ) 

4.  ,hen  the  complex  nefracli  on  index  of  the  liquid  is  calculated. 
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